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Abstract. - Here we study a pumping device capable of maintaining a density gradient and 
a flux of particles across a membrane. Its driving mechanism is based on the flashing ratchet 
effect powered by the random telegraph process in the presence of thermal fluctuations. Unlike 
Brownian motors, the concentrations at both reservoir boundaries need to be implemented as 
boundary conditions. The residence transition rates of the dichotomic flashing are related to the 
binding and hydrolysis of ATP molecules. The model is exactly solved and explored. The pump 
energetics is discussed and the relevant parameter values are tuned within a biological scale. 
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' Introduction. Active transport processes are ubiq- 
uitous in living matter. In fact, the cell cannot survive 
when specific out-of-cquilibrium tasks inside cease, mainly 
carried out by proteins such as kinesin, RNA-isomerase, 
dynein or rotatory ATP-ase. For instance, in order to 
rnaintain a suitable osmotically regulated state, the pump- 
ing of ions across the cell membrane is a crucial job that is 
performed by molecular pumps, which are, in turn, pow- 
ered by the hydrolysis of ATP molecules. Stemming from 
lihe studies on active transport in noisy environments [1], 
the stochastic features of pumps were originally addressed 
ip [2, 3] and applied to chemically driven electron pump- 
ing [4] . An experiment on artificial ion pores showing recti- 
fication was presented in [5] , where the observations where 
interpreted by a model containing flashing ratchets. Also 
ifi [6], experimental realizations of synthetic nanoscopic 
devices transporting potassium ions against their concen- 
tration gradient when stimulated with external field fluc- 
tuations were built. Recently it has been directly demon- 
strated that the ratchet mechanism is indeed of use in 
a biological system where molecular transport driven by 
fluctuations occurs in a membrane channel [7] . Other the- 
oretical works have studied the gating kinetics of ionic 
channels, based on the effect of non-equilibrium fluctua- 
tions combined with the ratchet effect [8]. A Brownian 
device with cyclic steps accounting for its conflgurational 
changes has been shown to reproduce the main features of 
experiments on the Na,K- ATP-ase ion pump [9]. Moving 
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to more simplistic descriptions, the study of active trans- 
port in an idealized Brownian pump focusing on the con- 
centration gradient created and maintained at both sides 
of a membrane has been analytically and numerically dis- 
cussed in [10]. 

Despite its relevance for modeling molecular pumps, the 
ratchet effect and its transport features against a concen- 
tration gradient have drawn very little attention. In this 
Letter we introduce a simple, yet revealing, two-state theo- 
retical model for a machine which pumps particles against 
a concentration gradient. The model is inspired by re- 
cent experimental works where the operating cycle of a 
channel has been shown to involve two steps (opening and 
closing) related to the ATP binding and hydrolysis respec- 
tively [11,12]. The pumping device we present consists of a 
ratchet potential, embedded in a membrane and bounded 
by particle reservoirs, which exhibits dichotomic fluctua- 
tions controlled by the ATP concentration. We do not 
focus on the precise structural details of the pump, but do 
it at the level of energetics. We consider the asymmetric 
ratchet profile as an energy barrier. In the steady state, 
we exactly calculate the concentration ratio between both 
reservoirs and the particle flux. We relate the performance 
of the device with ATP concentration, study its energetics 
and finally place it in a biological context. 

The model. Let us consider a pumping mechanism 
based on a flashing potential embedded in a membrane of 
length L, whose boundaries are infinite reservoirs of parti- 
cles with fixed mean densities po and pi . The dynamics of 
a particle in the membrane is determined by the following 
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Langevin equation in the over-damped regime, 



7i = -C/'(x,t)+e(i), 



(1) 



where U{x, t) is a time dependent potential (the prime de- 
notes position derivative) and ^(t) the random force aris- 
ing from the thermal fluctuations. It has zero mean and 
its autocorrelation is given by the fluctuation-dissipation 
theorem, 

{mm)^'ilkBT5{t-t'), (2) 

where 7 is the friction coefhcient and T the temperature of 
the environment. Whereas in the reservoirs the particles 
can diffuse freely, in the membrane they experience the 
potential 

U{x,t)^V{x)at). (3) 

where V(x), as depicted in figure 1, is a piecewise linear 
asymmetric ratchet potential: 



Vi{x)^Vo — , X G [0,(5L] 
oL 



V2{x) = 



L ^ X 
(1^' 



X e [5L,L\. 



(4) 
(5) 



Vq is the energy barrier height and 5 controls the spatial 
asymmetry by taking values between zero and unity. 

The flashing modulation C,{t) is a stochastic forcing in 
time which randomly switches the ratchet potential be- 
tween two states [1,2]. It is distributed according to a 
random telegraph (dichotomous Markov) process, flipping 
between values Ca = 1 (closed state) and — ^ (open 
state), with residence probability transitions Wa and Wb 
respectively [13,14]. Its mean value and correlation are 



(AC(t)AC(i')) 



Wa + Wb 

WgWb 

{Wa + WbY 



(6) 

-(Wa.+Wi,)\t-t'\ I'y'j 



where AC(i) = C,{t) — {({t)). While the asymmetric poten- 
tial breaks the spatial symmetry, the flashing transitions 
between states are responsible for the breaking of detailed 
balance in the system. 

Transition rates and ATP concentration. The 
hydrolysis of ATP molecules provides the energy input. 
It creates the necessary out-of-equilibrium conditions to 
drive molecular pumps. As a first level of approxima- 
tion to include the effect of ATP consumption into the 
model, we relate the ATP concentration, its binding and 
its hydrolysis processes to both residence transition rates. 
This two-state picture has been described in experiments 
in channels [11, 12] and it is a valid paradigm for the 
mathematical characterization of mechano-chemical en- 
ergy transduction. We consider that, when the pump is 
in the closed state (Ca = 1), it is the binding of an ATP 
molecule to a certain pocket center of the pump which in- 
duces conformational change to the open state {(b = 0). 
The transition rate Wa is thus dependent on the ATP 




Figure 1: Pump model: an asymmetric potential V{x) flashes 
its shape at random between a blocking (closed) and a flat 
(open) configuration with transitions rates Wa and Wb- 



concentration available in the surroundings, which can be 
modeled Michaelis-Menten kinetics law: 



wq- 



[ATP] 



Wq 



[ATP] 1 + 1/cr' 



(8) 



where kM is the so-called affinity constant and [ATP] is 
the concentration of ATP molecules. Only their ratio is 
relevant: a = [ATP]/kM. The rate Wb is given by the 
inverse of the machine's intrinsic time, which we consider 
constant: 

Wb = wq. (9) 

In the limit of a saturating [ATP], both transition prob- 
abilities are equal {wa = Wb = wq) indicating that the 
pump cannot operate faster. Medium [ATP] values imply 
a longer average time in the closed state. At low concen- 
trations. Wa is small and the pump seldom starts the cycle. 
Diffusive leak losses dominate and its pumping capacity is 
then very low. See scheme in figure 2. 

Numerical simulations. As a first test to show that 
the model just introduced can indeed create and maintain 
a concentration gradient, we use the simulation framework 
discussed in [10]. This method allows for the numerical im- 
plementation of the reservoir boundary conditions, which 
is not trivial. We then simulate non-interacting particles 
following the Langevin equation (1), under thermal white 
noise and the random telegraph process. In figure 3, we 
plot the histogram of the position of the particles, be- 
ing illustrative of the density profiles typically obtained 
in the steady state at zero net flux. It is clear that the 
pump is able to generate and maintain a concentration ra- 
tio (pi po), while the edges emulate particle reservoirs 
with uniform concentrations. 

Theoretical analysis. Consistent with the Langevin 
equation (1), the corresponding partial differential equa- 
tions for the time evolution of the concentration of parti- 
cles in each flashing state. Pa = Pa{x, t) and Pb = Pb{x,t), 
originate from a combination of the Fokker-Planck de- 
scription for the continuous variable x and the Master 
equation for the discrete states a and b. They read 
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Figure 2: Illustrative representation of the effect of ATP con- 
centration in the switching rates of the dichotomic flashing 
ratchet mechanism. 




Figure 3; Density profile obtained directly from numerical sim- 
ulations of the Langevin equation in the steady state at zero 
flux. 



[2,13,14] 

dtPa = J-^d:,[V'{x)+kBTd,]Pa-WaPa+WbPb, 
dtPb = J-^d^[kBTd^]Pb~WbPb + WaPa. (10) 

It is convenient to work with dimensionless quantities. 
Wc first re-scale space and time: z = x/L and s = 
t{kBT)/{-fL'^). Consequently, the new transition rates are 
/ = w{'jL'^)/ (ksT) (for Wa, Wb and wq) and L is absorbed 
in the concentration distributions. We introduce the fol- 
lowing densities: the total density of particles p = Pa + Pb 
and the auxiliary function n = Pa — Pb- In the steady 
state, eqs. (10) are governed by a new set of equations for 
each region (indicated by the subindex i = 1,2): 



-J, (11) 
4/fc + /a)-p(/f>-/a), (12) 



where Fi = -vo/2S and F2 = vo/2{l - 6). vq = VolksT 
measures the relative strength of the potential with re- 
spect to the thermal energy. J is the physical total flux of 
particles across the membrane. 

Substituting (11) in (12) one gets a linear ODE for Pi{z), 

p;'{z)-2F,p;{z)-{fa + fb)pAz) + 2fbF,p,{z) = J{fa + h), 

(13) 

whose solution has the standard form 

p,{z) = Cie^-^ + C,2e^-^ + ase^'^' + {J/F,)a, (14) 

where a = {fa + fb)/2fb- The constants Cy are, so far, 
unknown parameters. Instead, the constant coefficients 
Xij are found by inserting pi{z) back in (13) and solving 
the algebraic equation 

\% - 2F,\l - {fa + fb)Kj + 2fbF, = 0. (15) 

Plugging (14) in equation (11) and defining Sy = {Xij / Fi~ 
1) and (3 = {fb - fa)/2fb, one finds 

7r,(z) = saC.ie^'i" + s^sC.se^'^" + s^sase^'^' + {J/F^)(3. 

(16) 



Boundary conditions. The formal solution of the 
concentration of particles in the steady state depends on 
seven unknown constants (the coefficients C.y and the flux 
J), whose value is determined by boundary conditions. In 
contrast with models for Brownian motors, in the present 
Brownian pump one should not impose normalization nor 
the periodicity condition for the probability. Instead, 
the total concentrations at left and right boundaries, po 
and pi respectively, are externally fixed: pi(0) = po and 
P2(l) = pi- At the boundaries (in contact with the par- 
ticle reservoirs) the state probabilities can be factorized 
as the total concentration times the stationary weight 
the state, complying with Pa{0) Pofb/{fb + fa) and 
Pb{0) ^ Pofa/{fb+fa) (and analogously at z = 1). Rewrit- 
ten in terms of 7r(z), this supplies two new conditions: 
7ri(0) = rpo and 772(1) =rpi, where r= {fb- fa)/ {fb+ fa)- 
Third, p{z) and 7r(z) must be continuous functions con- 
necting both regions: pi{5) = P2{S) and Tri{S) = 7r2(<5). 
Last, the total flux J associated to p{z) has implicitly been 
assumed continuous across both zones. The same holds for 
the one associated to Tr{z) in (12), which in practice yields 
to impose continuity on tt[{z) — Fi[pi{z) + ni (z)] . This last 
boundary condition can be recast as a property of the 7r(z) 
derivatives: (F2 - Fi)[pi{S) + tti{S)] = 7r^((5) - tt[{5)- 

Solution. The above cumbersome boundary conditions 
can be compactly written as a system of linear equations: 

see eq.(17), 

from which the flux J (and also the density profile p{z) by 
means of Cij ) can be obtained as a function of the system 
parameters. New constants arc defined for compact nota- 
tion: AT = {I/F1-I/F2) and % = K, (s^ + F2/F1 - 1). 
In the zero flux case (J = 0) such linear system has to be 
rearranged. The ratio of concentrations pi/po, which is 
then the relevant observable, is easily found by imposing 
a minor of the coefficient matrix equal to zero (ensuring 
that the solution of the problem exists): 

see eq.(18). 
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Parameter exploration: zero flux case. When the 
pump is forced to operate at vanishing particle flux, pi 
and Po become boundary conditions, whose ratio is the 
relevant quantity to explore as a function of the parame- 
ters of the system. When examining the behavior of pi/po 
as a function of the asymmetry parameter S, we check that 
for (5 = 0.5 (spatially symmetric system) no net pumping 
is possible, thus leading to pi = po- Moreover, the solu- 
tion of the problem is anti-symmetric when i5 and 1 — d 
arc interchanged. The pumping ability clearly depends 
on the relative energetic barrier vq- For low values of vq, 
the ratchet can barely pump the particles because diffu- 
sive losses dominate, whereas as the barrier height is in- 
creased, the concentration ratio does so. For very high 
values of vq, the fluctuating potential does not improve its 
effect anymore and its pumping capacity saturates. All 
these features were verified (data not shown) and used as 
consistency tests for the analytical solution. 

The ratio pi/po as a function of the transition rates of 
the dichotomic noise is explored in figure 4. By means 
of varying /o we observe that there is an optimum value 
which maximizes the particle gradient. For a very fast 
switching (/o high), the particles do not have time to dif- 
fuse in the open configuration and cannot be biased in 
the closed configuration. The opposite limit is also phys- 
ically intuitive. In the adiabatic case (/o small), no par- 
ticle gradient can be created: the flashing time-scale is 
much greater than diffusion, which prevents any possible 
separation of particles induced by the ratchet effect. Fur- 
ther adjustments on 6, vq and a may allow to obtain even 
higher concentration ratios, optimizing the performance of 
the pump. Figure 4 depicts how tuning both parameters 
fa and fb (equivalently we vary /o and a) can increase the 
ratio pi/po- For very small [ATP], fa tends to zero. The 
ratchet is in the closed state and it takes a very long to 
open. This leads to pi po. For high values of a, the 
rates are equal, fa = fb = fo, and /o alone controls the 
pumping capacity. In the intermediate region, there is a 
maximum distinctive of an optimal behavior, since a can 
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Figure 4: Concentration ratio pi/po versus the characteristic 
frequency /o for different a values [5 = 0.33 and vo = 10). 

independently improve the pumping capacity given a fixed 
value of /q. 

Parameter exploration: non-zero flux case. In 

figure 5-top we plot the flux J as the concentration pi is 
increased (given a fixed po = 1) for several values of vq. 
A restoring entropic force appears opposing the pumping 
and leading to a decrease of the total fiux. The decrease 
of the fiux is linear. When the barrier tjq is higher, the 
decay is softer due to a more powerful functioning of the 
pump (at the expense of more energy) . Even if the pump 
is able to carry a non-zero net flux of particles against the 
concentration gradient, if pi is sufficiently increased, the 
particle flux can be reversed. 

The analytical prediction of the flux J as a function of 
the rate /o for several values of pi (where po = 1) is shown 
in figure 5-bottom. When pi = po and /o = 0, the device 
is stopped and the net flux vanishes. As the pumping rate 
is increased, transport starts. The higher pi, the lower 
the flux J, even becoming negative for great concentration 
values. In that case, if the flashing rate is increased, the 
machine can reverse the sign of the flux back to J > so 
that particles are then pumped against the concentrations 
imposed by the reservoirs. Finally, for very high rates the 
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Figure 5: Top: Total flux J as a function of concentration p\ 
when the effective height of the barrier vo is varied (po = 1, 
5 — 0.33, fo = 100 and a = 0.5). Bottom: Flux J as a function 
of the flashing rate fo at different concentrations pi (po — 1, 
S = 0.33, vo = 20 and a = 0.5). 



flux steadily tends to zero because, when the flashing is 
almost instantaneous, there is no effective pumping. 

Energetic characterization. The energetic aspects 
are relevant and abundant in the literature of motor mod- 
eling [15-17] and should definitely be addressed in their 
pump counterparts. A suitable definition of efficiency for 
pumps is not obvious. It is unclear how to interpret and 
quantify the payoff achieved from the pumping process as 
a whole. However, from a purely energetic point of view, 
meaningful quantities can still be defined consistently. 

First, the input source of energy can be conveniently 
studied as follows. Working in the overdamped approxi- 
mation, every time the potential is lifted (due to the 
modulation) energy is injected into the system. This en- 
ergy is used to actively create the concentration profile 
and maintain it. When the pump fiashes down to the 
open configuration, no energy is injected nor recovered. 
In such situation, the pumping energy input per unit of 
time is mathematically quantified as 

Ein^ I dxV{x)wbPb{x). (19) 
Jo 

The mean power input that the dichotomic flashing in- 
serts into the system is a controllable feature of the model 
by a judicous choice of its pameter values. An explicit 
expression for the above quantity (in dimcnsionlcss vari- 



ables) can be obtained by splitting it into two parts, each 
corresponding to the piecewise linear potential: 

^in _ (' [pi(z)-^i(z)]/2 , r\ [p2{z)~^2{z)]/2 

vofo^Jo S/z (l-5)/(l-z) ■ 

(20) 

By using the state probabilities (14) and (16) we are left 
with sums of exponentials, whose integrals are trivial. In 
flgure 6 we show the behavior of Ei^ at zero flux as the 
typical frequency /o is varied (solid line) for several values 
of vq. The faster the flashing, the more energy per unit of 
time is inserted in the system. Similarly, when the barrier 
Wo is increased, the power input rises very quickly. 

Secondly, the power output produced by the machine 
can be quantifled as the product i?chom = J where J is 
the particle current and A0 = ksThi^pi/pQ) is the chemi- 
cal potential difference between both reservoirs. Note that 
at zero flux, the power output vanishes, namely, the pump 
is using energy but it is not performing work. This is 
equivalent to the stalling regime in motors, where the en- 
ergetic efficiency is zero since there is consumption but no 
work performance. 

On the whole, the benefit pursued by a pump is often 
not an energetic one. For instance, maintaining a con- 
centration gradient at zero flux can indeed be very im- 
portant for biological functions. The energetic input can 
still be used to constrain the biophysical tuning of the 
model, which we address as follows. In connection to real 
biomachines, one ought to require the input energy per 
cycle not to be greater than the energy released from the 
hydrolysis of one ATP molecule, which is approximately 
20fcsT'. Given that the typical period of the pump cycle 
is Tcycic = + w^^ , the power released by the ATP is 
Eatp = -E'ATp/''cycic which, in dimensionless units, yields 

20 

Eatp = fo^r—^ > i?i„. (21) 

2 + l/CT 

In flgure 6 we plot the above expression (dashed line) as 
a function of the frequency /q. 

Biophysical scale. So far we have solved the equations 
of the model and examined the behavior of the pump as 
a function of dimensionless parameters. Can all that be 
placed in a biophysical scale [18-21]? We use a tuning 
approach that combines constraints and optimization cri- 
teria. Let us consider the length of the molecular pump 
of the order of i ~ 5nm and an effective opening section 
of 5 ~ Inm^. At biological conditions (T = 310K) the 
thermal energy is fcsT ~ 4.28 • 10~^^J. Given an effective 
viscosity inside the membrane two orders of magnitude 
bigger than that of the water, the resulting diffusion con- 
stant of the medium yields D = ksT/^ ~ 10~^^to^/s. 
This leads to a typical time-scale associated with diffu- 
sion of the order r — LF' /D ~ lO^^s. Then, when the 
optimal fluctuating rate /o found in our explorations is 
mapped to dimensional units we flnd wq ~ 10^ ^ lO'^Hz, 
which is compatible with typical reported rates of func- 
tioning of the Na,K-ATP-ase pump under the effect of 
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Figure 6: Input power Ein as a function of tiie flasiiing rate /o 
for three different values of vq {S = 0.33, a = 0.5, po = 1 and 
J = 0). The dashed hne represents the power related to the 
hydrolysis of one ATP molecule per cycle. The parameter space 
of the model is tuned at the black dot according to biophysical 
constraints. 



an oscillating electric field [22]. Guided by the energetic 
calculations developed previously, we impose that the en- 
ergy spent per cycle is close but not higher than the one 
obtained from the ATP hydrolysis. See the black dot in 
figure (6). This sets the possible space of parameters. For 
instance, a plausible potential barrier Vb is expected to be 
of the order of a few tens of fcsT. Higher values can lead to 
greater concentration ratios but their energetics is not bio- 
logically meaningful. In this way, by recursively adjusting 
the parameter space where the pump functions optimally, 
bearing in mind biological constrains, the model is shown 
to be consistently translated into dimensional units of the 
order of those of a real molecular pumps. 

Final remarks and conclusions. The interplay be- 
tween experimental observations and theoretical models 
substantially helped to understand the basic features of 
molecular motors. We believe that applying now the basic 
mechanisms of the ratchet effect to molecular pumps will 
contribute to understand the fundamental working prin- 
ciples of protein machines embedded in cell membranes. 
Although the model presented here is a dramatic sim- 
plification of the great complexity and richness of a bi- 
ological pump, it contains many interesting and new el- 
ements that capture its basic biophysical aspects. First, 
the explicit implementation of concentrations as bound- 
ary conditions at both ends of the membrane, leading to 
the characteristic quantity pi/po- Second, the use of the 
simulation framework developed in [10], which is prone 
to be applied in models whose an analytical solution is 
not possible. Third, a simple coarse-grained connection 
between transition rates and ATP binding and hydrolysis 
processes. Fourth, the energetic characterization of the 
model and its connection to biophysical values. Finally, 
it is worth commenting on the straightforward extension 
of this approach for particles of charge g in a membrane 
potential $ by including an electrostatic force Fq = /L 



in the equations of motion. The pump then needs to sup- 
ply extra energy, reducing the flux and the concentration 
ratio accordingly. 

On the whole, we expect that the present work will be 
of use as a starting point to make progress, from a phys- 
ical point of view, towards future investigations of active 
processes in membranes. 
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